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POSSIBLE  ROLE  OF  VIBRATIONALLY  EXCITED  HYDROGEN  IN  THE 
CHEMISTRY  OF  THE  EARTH’S  UPPER  ATMOSPHERE 

I.  Introduction 

The  chemistry  of  vibrationally  excited  molecules  in  the  earth's 
atmosphere  is  a  subject  which  has  received  substantial  attention  (Dalgarno, 
1963;  Cvetanovic,  1974;  Taylor,  1974,  Vlasov,  1976).  Much  of  this  attention 
has  been  due  to  the  fact  that  emission  from  vibrationally  exlcted  molecules  in 
the  upper  atmosphere  is  known  to  be  responsible  for  major  portions  of  the 
earth's  airglow,  most  notably  the  Meinel  bands  of  OH  (Bates,  1982  and 
references  therein).  Infrared  emission  from  vibrationally  exicted  molecules 
in  the  thermosphere  is  also  expected  to  play  a  major  role  in  atmospheric 
cooling  (Gordlets  at  al. ,  1982). 

There  is  little  evidence  that  chemical  reactions  of  vibrationally  exlcted 
species  play  an  Important  role  in  controlling  the  net  photochemistry  of  the 
earth's  upper  atmosphere,  however.  By  this,  we  mean  that  it  is  not  expected 
that  processes  of  the  type 

A*  +  B  ♦  C  +  D  (la) 

will  be  important  while  processes  such  as 

A  +  B  -»•  C  +  D  (lb) 

are  not  (an  asterisk  is  used  here  to  represent  vibrational  excitation). 
Crutzen  (1974)  suggested  one  case  where  vibrational  excitation  could  play  a 
very  Important  role.  This  case  is  the  reaction 

0(^P)  +  >  OH  +  H  (2a) 

which  is  known  to  occur  very  slowly  at  low  and  moderate  temperatures  due  to  a 
large  (  8.5  kcal/mole)  activation  energy.  Because  the  H2  vibrational 
quantum  is  large  (10.2  kcal/mole)  and  vibrational  excitation  of  H2  is  likely 
in  H2  forming  reactions,  Crutzen  suggested  that  reaction  2  could  be 
sufficiently  accelerated  by  vibrational  excitation  of  H2  that  it  might  become 

\fanuscript  approved  September  7,  1983. 
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one  of  Che  dominant  loss  processes  for  H2*  He  showed  that  If  one  assumed  a 

—  1  1  —  1 

temperature  independent  rate  constant  of  2.9  x  10  cm  molec  for  the 
reaction 


0(^P)  +  H^Cv  >1)  +  OH  +  H  (2b) 

the  predicted  concentration  of  H2  at  and  above  the  80  km  level  of  the 

atmosphere  decreased  by  nearly  50%. 

This  increased  destruction  of  H2  by  reaction  (2b)  led  to  a  substantial 
increase  in  the  mixing  ratios  of  H2O  above  70  km  calculated  by  Crutzen 
(1974).  Odd  hydrogen  (HO^)  formed  as  a  result  of  H2O  photolysis  is  known  to 
lead  to  ozone  (O3)  depletion  by  catalytic  cycles,  and  such  reductions  are 

expected  to  be  especially  strong  at  higher  altitudes  (Prather,  1981).  As  the 

ability  to  measure  0^  in  the  upper  mesosphere  and  lower  thermosphere  is 
increasing,  the  uncertainty  in  the  H2O  concentration  of  the  magnitude 

it 

indicated  by  Crutzen  (1974)  due  to  the  possible  role  of  becomes 

unacceptably  large.  Thus,  it  is  important  that  the  role  of  reaction  (2b)  as  a 
means  of  altering  the  upper  atmospheric  H2  and  H2O  balance  and,  therefore,  the 
O2  concentration  be  established 

Recent  experimental  work  (Light,  1978)  indicates  that  Crutzen  seriously 
overestimated  the  effect  of  H2  vibrational  excitation  on  the  rate  of  reaction 
2,  and  theoretical  calculations  (Schinke  and  Lester,  1979)  suggest  that  more 
than  two  vibrational  quanta  of  excitation  are  needed  for  the  rate  of  reaction 
2  to  even  begin  to  approach  the  value  suggested  by  Crutzen. 

In  this  report  we  examine  in  some  detail  the  aeronomy  of  vibrationally 

* 

excited  H2,  hereafter  referred  to  as  H-  ,  in  the  earth's  upper  atmosphere.  We 

consider  processes  responsible  for  production  and  destruction  of  H2  ,  and 

attempt  to  make  reasonable  estimates  for  the  vibrational  state  specific 

chemical  reaction  and  vibrational  quenching  rates.  Finally,  we  combine  the 

expressions  obtained  with  a  knowledge  of  the  photochemistry  of  H2  in  the 

atmosphere  (Hunten  and  Strobel,  1974;  Liu  and  Donhue,  1974;  Schmidt,  1974; 

* 

Koshelev,  1976;  Penner  et  al. ,  1977)  to  conclude  that  consideration  of  H2  is 
not  expected  to  appreciably  alter  the  calculated  upper  atmospheric  H2 
densities. 
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II.  Chemical  Model 


In  this  section  we  will  briefly  review  all  chemical  reactions  which  we 

* 

expect  to  be  important  in  the  aeronomy  of  ;  we  will  defer  estimation  of 

state  specific  rate  constants  to  Section  III.  We  will  consider  processes 

responsible  for  production,  loss  by  reaction,  and  loss  by  inelastic  processes 

* 

of  • 

A.  Production  of  H* 

Crutzen  (1974)  suggested  three  processes  which  can  reasonably  be 

A 

assumed  to  give  rise  to  H2  in  the  earth's  upper  atmosphere: 

H  +  HO2  -  H2(v)  +  O2  (3a) 

H2O  +  hu  >  H2(v)  +  0(^D)  (4a) 

H  +  H  +  M  +  H2(v)  +  M  (5) 

Of  these,  the  small  concentrations  of  H  in  the  earth's  upper  atmosphere  and 
the  slow  rate  for  R5  combine  to  make  it  unimportant.  Model  calculations  of 
Hunten  and  Strobel  (1974)  show  that  the  fractional  contribution  of  R5  to 
H2  production  near  its  peak  at  80  km  must  be  not  greater  than  1  part  in  10^. 
Thus,  we  will  restrict  our  attention  to  the  potentially  important  reactions  R3 
and  R4a. 

We  note  that  R3a  and  R4a  represent  each  only  one  possible  product 
path  for  the  respective  reactions.  Alternate  (and  frequently  more  prevalent 
paths)  are 


li  +  HO2  ->•  2CH  (3b) 

H2O  +  0  (3c) 

H2O  +  hu  -*■  H  +  OH  (4b) 

-*■  0  +  H  +  H  (4c) 


While  the  exact  portion  of  reactions  3  and  4  resulting  in  H2  formation  by 
pathways  3a  and  4a,  respectively,  is  not  precisely  determined  (Hack  £t_  al . 
1978;  Sridharan  et  al.  1982;  Stief  et  al. ,  1975;  Slanger  and  Black,  1982);  the 
results  of  our  analysis  do  not  depend  strongly  on  the  branching  ratios  R3a/R3 
and  R4a/R4  unless  we  have  greatly  underestimated  them. 
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B.  Loss  of  H*  by  Chemical  Reaction 

Three  H2  destruction  reactions  were  considered  by  Crutzen  (1974): 


0(^P)  +  H2 

♦  OH  +  H 

(2) 

0(^D)  +  H2 

->■  OH  +  H 

(5) 

OH  +  H, 

+  H2O  +  H 

(6) 

Of  these,  it  is  expected  the  R2  is  the  only  one  which  will  be 

•k 

important  for  .  Reaction  5  is  known  to  occur  rapidly 
(k^  ■  1.0  X  10”^  cm^molec'^sec""^  and  no  temperature  dependence) 
with  ground  vibrational  state  H2  1982),  so  vibrational  excitation  is  not 

expected  to  appreciably  accelerate  it.  Supporting  this  expectation,  quasi- 
classical  trajectory  calculations  for  R5  show  only  a  10%  increase  in  reaction 

cross  section  on  H2  vibrational  excitation  (Schinke  and  Lester,  1980),. 

^  1  .  ,  3  ,  * 

Since  [0(  D)1  <<  ]  0(  P)1  ,  we  may  thus  safely  neglect  the  effect  of  R5  on  II2 

Reaction  6,  which  is  known  to  have  a  substantial  (4  kcal/mole) 

activation  energy  (JPL,  1982)  could  be  expected  to  be  appreciably  accelerated 

by  H2  vibrational  excitation,  especially  at  low  temperatures.  At  room 

temperature,  single  quantum  excitation  of  H2  is  known  to  increase  the  rate  of 

R6  by  roughly  a  factor  of  100  (Glass  and  Chaturvedi,  1981),  corresponding  to 

an  approximate  decrease  in  the  activation  energy  from  4  to  2.5  kcal/mole 

(Truhlar  and  Isaacson,  1982)  Even  if  multiple-quantum  excitation  of  H2 

sufficiently  accelerates  R6  that  there  is  no  effective  activation  barrier,  the 

low  concentration  of  OH  in  the  mesosphere,  as  opposed  to  that  of  0(  P) , 

(Prather,  1981)  indicates  that  R2  should  be  much  more  important  than  R6,  and 

A 

we  will  thus  consider  R2  as  the  only  reactive  loss  process  for  H2  . 

C.  Loss  of  H*  in  Inelastic  Collisions 

We  will  need  to  consider  three  major  non-reactive  loss  process 

for  H2  : 


H2(v)  +  N2  ^  H2(v'  <  v)  +  N2  (7) 

H2(v)  +  02-*-  H2(v  <  v)  +  O2  (8) 
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H2(v)  +  0 


(9) 


+  H^Cv  <  v)  +  0 

■k 

where  we  write  here  as  H^Cv)  to  reinforce  the  state-specificity  of  the 
quenching  rates.  Since  there  is  no  information  on  the  rate  of  R8,  we  will 
assume  that  the  rates  for  quenching  of  H2(v)  by  ^^2  and  0^  are  identical.  This 
is  a  reasonable  assumption,  given  the  similarities  of  N2  and  O2  in  terms  of 
mass,  vibrational  frequencies,  and  Lennard- Jones  potential  parameters 
(Lambert,  1977). 

The  only  measurement  of  the  rate  of  R7  is  for  v>l  at  room  temperature 
(Bott,  1976)  and  is,  therefore,  not  directly  applicable  to  the  high  excitation 
and  low  temperatures  encountered  in  the  upper  atmosphere.  Ue  will  attempt  to 
make  reasonable  approximations  for  the  v  and  T  dependence  of  the  rate  of  R7. 
There  are  no  measurements  of  the  rate  of  R9,  and  no  results  were  presented  for 
H2  vibrational  deactivation  in  quasi-classical  trajectory  calculations  on  R2 
(Schinke  and  Lester,  1979),  so  we  will  be  forced  to  estimate  the  rate  of  R9. 


III.  Estimation  of  State-Specific  Reaction  Rates  and  Branching  Ratios 
A.  Production  of  H* 

1.  H  +  HO,  -  H,  (v)  +  0, 

d  A'  li 

Reaction  3a  is  sufficiently  exothermic  (2.25  eV  if  the  O2  is 

produced  in  its  ground  electronic  state)  that  the  v  »  0-4  vibrational  levels 

of  H2  are  energetically  accessible.  If  the  O2  is  produced  in  one  of  its  low- 

lying  electronically  excited  states  (a^A  ,  b^Z^)  ,less  energy  is  available  for 

S  S 

H2  product  internal  excitation.  Experiments  indicate  that  the  yield  of 
electronically  excited  O2  in  reaction  3a  is  quite  small  (Washida  et  al. ,  1978) 
and  we  will  assume,  therefore,  that  reaction  3a  proceeds  only  to  ground 
electronic  state  products. 


No  information  is  available  on  the  distribution  of  internal  energy  in  the 
products  of  reaction  3a.  Quasi-classical  trajectory  calculations  on  a  related 
reaction 


H  +  HCO  +  H2  +  CO  (10) 
have  been  performed  (Fanantos  and  Murrell,  1980).  They  showed  that  between 
25%  and  30%  of  the  available  product  energy  goes  into  H2  vibration  in  the 
temperature  range  300-1000K,  and  that  the  H2  vibrational  state  distribution  is 
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approximately  Boltzmann.  If  we  assume  these  conclusions  carry  over  to 
reaction  3a,  we  estimate  a  product  state  distribution  for  v  «  0,  I,  2,  3,  and 
4  respectively  of  0.35:  0.25:  0.18:  0.13:  0.09. 

2.  H7O  +  hv  *  H7(v)  +  0(^D) 

Nothing  Is  known  of  the  H2  vibrational  state  distribution  in  this 
process.  Substantial  vibrational  excitation  of  H2  is  energetically  permitted; 
the  endothermic! ty  of  the  dissociation  process  is  7.06eV,  and  a  photon  at 
1216  (Ly-a  ,  expected  to  be  the  most  Important  vacuum  ultraviolet 
wavelength)  contains  10.20eV,  leaving  3.14eV  available  for  product 
excitation.  Given  that  there  is  an  appreciable  change  in  the  H-H  distance  on 
photodissociation,  substantial  vibrational  excitation  of  the  H2  product  is  a 
distinct  possibility.  If  we  estimate  that  20Z  of  this  available  energy  goes 
into  H2  and  assume  that  the  vibrational  states  are  populated  with  a  Boltzmann 
distribution,  a  product  state  distribution  for  v  ■  0,  1,  2,  3,  4,  and  v  >  5  of 
0.4:  0.2:  0.1:  0.08:  0.05:  0.06  is  expected.  This  distribution  is  reasonably 
close,  especially  at  the  high  vibrational  excitation  end,  to  that  which  would 
be  expected  on  a  purely  statistical  basis  (i.e.  product  states  are  populated 
proportionately  to  the  density  of  states). 

B.  Loss  of  H*  by  0  (3?)  +  H^Cv)  >  OH  4-  H 

The  available  data,  both  experimental  (Light,  1978)  and  theoretical 

(Schinke  and  Laster,  1979)  indicate  that  H2  vibrational  excitation  leads  to  an 

enhancement  in  its  rate  by  a  factor  much  smaller  than  that  assumed  by  Crutzen 

(1974).  At  room  temperature  Light  (1978)  measured  a  rate  constant  of  1.0  x 

10~^^  cm^molec”^sec”^  for  R2b  with  v=l.  No  measurements  of  this  reaction  at 

the  lower  temperatures  of  the  earth's  mesosphere  and  lower  thermosphere  have 

been  reported.  Two  quantum  excitation  of  H2  has  been  calculated  to  lead  to  an 

additional  factor  of  100  enhancement  over  the  rate  for  R2b(v  ■  1)  at  300K 

—  1  9 

(Schinke  and  Lester,  1979),  which  gives  an  approximate  value  of  10 

"i  _i  _1 

cm  molec  sec  for  R2b  (v  *  2)  at  300K.  In  order  for  R2b  to  occur  rapidly 
(k^^^  >  10-1 ^cm^molec  ^  ^e  clearly  require  high  vibrational  excitation. 

As  we  showed  earlier,  population  to  v  »  5  is  probable,  so  detailed 
consideration  of  these  higher  levels  is  appropriate. 
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A  simple  picture  which  permits  one  to  make  reasonable  predictions  for 
the  effect  of  vibrational  energy  on  chemical  reactions  involves  a  combination 
of  transition  state  theory  and  the  vibrationally  adiabatic  theory  of  chemical 
reactions  (Kuppermann,  1979  and  references  therein).  In  this  picture,  one 
specifically  accounts  for  the  fact  that  large  effective  energy  barriers  to 
chemical  reactions  may  persist  even  when  the  amount  of  vibrational  excitation 
is  very  large  (Kaye,  1982).  This  method  has  been  applied  to  R2  (v  ■  0,1)  and 
shown  to  give  reasonable  results  (Lee  et  al. ,  1982b),  and  we  use  a  simplified 
version  of  this  model  to  make  predictions  for  the  temperature  and  quantum 
state  dependence  of  R2b. 

The  crucial  step  in  this  method  is  the  recognition  that  the  position 
and  energy  of  the  transition  states  (represented  by  t  )  of  the  chemical 
reaction 


0  +  H2(v)  -  O.H.H*(v)  0H(v)  +  H  (11) 

depend  on  the  reagent  vibrational  quantum  number  v.  The  transition'-state 
theory  expression  for  the  state-specific  is 


k  (T) 

V 


K(T) 


kT 


tr 


£?  EfZ  (E  ,)  "z 

tr  tr  '  rot  V 


exp(-aV  VkT) 


(12) 


2  3/2 

where  f^^  is  a  translational  partition  function,! »  (2TTmkT/h  )  ,  f^  is 
the  bending  partition  function  not  including  zero-point 
energy,  f^  »  il  ~  exp  (“hv^  ^  >  ^”^1  ^rot  ^  rotational  partition 
function.  is  the  v-dependent  height  of  the  vibrationally  adiabatic 
barrier  to  chemical  reaction  which  includes  stretching  and  bending  zero-point 
energies.  K(T)  is  a  dynamical  correction  factor  which  we  will  take  to  be 
unity  for  reasons  described  below.  In  usual  applications  of  transition  state 
theory  there  are  also  terms  for  reagent  vibration  and  complex  "symmetric- 
stretching"  type  modes  but  these  are  not  considered  here  because  our  model 
restricts  us  to  one  vibrational  state  at  a  time;  in  any  case  both  these  modes 
are  at  a  sufficiently  high  frequency  that  at  mesospheric  temperatures  their 
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value  is  unity. 

The  position  and  energy  of  the  vlbrationally  adiabatic  barriers  for 

R2  are  determined  by  the  method  dea.rlbed  by  Lee  et  al.  (1982)  except  that 

eigen-values  were  determined  nur  =‘rlcally  (Truhlar  1972)  rather  than  by  a  fit 

to  a  Morse  oscillator,  using  the  MOD  POL-CI  XMOS  potential  energy  surface  (Lee 

et  al.  ,  1982a).  Transition  state  rotational  partition  functions  were 

calculated  using  the  geometries  calculated  for  each  transition  state,  under 

the  rigid-rotor  approximation.  Ho  rotational  partition  functions  are 

independent  of  v  in  the  harmonic  oscillator  approximation  for  Ho.  Complex 

bending  partition  function  were  calculated  using  an  assumed 

form  ^  ■  A  exp(-bs^)l  for  the  dependence  of  the  bending  frequency  on  the 

distance  s^  of  the  v^^  transition  state  along  the  reaction  coordinate  from  the 

saddle  point.  Since  the  MOD  POL-  Cl  RMOS  surface  (Lee  et  al.,  1982a)  is 

defined  only  for  collinear  configurations,  we  estimated  b  using  the  Johnson- 

Winter  (JW)  potential  energy  surface  (Johnson  and  Winter,  1977).  A  was 

assumed  to  have  the  value  of  514  cm“^  determined  by  Walch  et  al.  (1980)  for 

■1 

their  A”  surface.  Bending  frequencies  were  calculated  using  standard 
formulas  for  linear  triatomic  molecules  (Johnston,  1966).  Values  of  various 
quantities  in  Eq.  12  and  parameters  used  to  calculate  these  values  in  Eq.  12 
are  shown  in  Table  1. 

We  note  that  the  MOD  POL-CI  RMOS  surface  does  not  have  the  correct 


asymptotic  vibrational  energy  levels  for  H2»  but  this  is  not  expected  to  lead 
to  appreciable  errors  except  for  v-0  (where  the  barrier  is  assumed  to  be  at 
the  saddle  point).  The  vlbrationally  adiabatic  barriers  are  sufficiently  far 
in  the  entrance  channel  (closer  to  0  +  H2  than  OH  +  H)  that  the  heights  and 
positions  of  the  barriers  are  probably  correct;  only  their  absolute  energies 
are  wrong,  and  they  do  not  enter  into  these  calculations.  The  dynamical 
correction  factors  K(T)  are  taken  to  be  unity  as  a  reasonable  compromise 
between  two  competing  effects.  At  low  T,  quantum  mechanical  tunneling  tends 
to  make  K(T)  large  (>1)  while  at  high  T,  reaction  probabilities  are 
always  <1,  so  K(T)  <  1  also.  Thus  we  feel  our  choice  for  K(T)  is  a  physically 
reasonable  one.  If  anything,  this  may  be  a  slight  overestimate  of  K(T),  which 
means  we  are  overestimating  k^(T)  and  thus  overestimating  the  importance  of 
R2b.  More  accurate  determination  of  K(T)  is  possible  (Bowman  et  al. ,  1981; 
Lee  al. ,  1982b)  but  becomes  expensive  because  of  the  large  Internal 
excitation  we  must  consider  here. 
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Table  1 


Values  of  Quantities  Used  In  Estimating  Rate  Constants  for 


the  Reaction  0  + 

H2(v)  X  oh  + 

H 

a 

b 

c 

d 

d 

a 

V 

X 

V' 

V 

v';h2(v)^^ 

X 

I” 

V 

± 

s,v 

X 

b,v 

X 

AV” 

V 

0 

-3.486 

-3.872 

2.33 

1.74 

3.22(-39) 

1636 

514 

0.450 

1 

-3.190 

-3.375 

3.01 

1.49 

4.36(-39) 

3750 

488 

0.246 

*0 

-2.812 

-2.913 

3.39 

1.46 

5.22(-39) 

4090 

451 

0.157 

3 

-2.429 

-2.485 

3.83 

1.44 

6.32(-39) 

4210 

389 

0.104 

4 

-2.048 

-2.091 

4.05 

1.43 

6.93(-39) 

4280 

355 

0.087 

5 

-1.703 

-1.732 

4.27 

1.44 

7.57(-39) 

4280 

319 

0.069 

a)  Potentials  are  measured  (In  eV)  with  respect  to  a  zero  energy  at  three 
infinitely  separated  atoms.  V'  ■  height  at  top  of  vibrationally 

adiabatic  barrier. 

VfH^Cv)]  *  energy  of  vibrational 
level  of  H2 

AV*  *  vibrationally  adiabatic  barrier 
height 


b) 

c) 


R”  is  xy  internuclear  distance  (in  bohr)  at  v*’"  transition 

9  t  h 

is  moment  of  inertia  (in  g  cm'^)  of  v*^"  transition  state. 
Numbers  in  parentheses  represent  powers  of  ten  by  which  the 


state 

non-enclosed 


d) 


numbers  should  be  multiplied, 
j*  .  ,  vj  „  are  "symmetric  stretch 


3  .V 


b,v 


respectively  (in  cm”^)  at  v 


.th 


and  bending  vibrational  frequencies 
transition  state 
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Using  Che  data  In  Table  1  and  eq.  12  we  find  an  expression  for  Che 


race  consCanC  of 

k^(T)  -  a^[l-exp(-byT)T^T^/^exp(-cyT) 


(13) 


Values  of  Che  consCancs  a^,  b^,  and  c^,  along  wlch  ky  aC  T-200,  300K  are  given 
In  Table  2.  The  300K  resulc  for  v  »  1  agrees  qulce  well  wlch  Che  resale  of 
l(+0.9,  0.6)  X  inolec~^sec~^  deCermlned  by  LlghC  (1978). 


C.  Loss  of  H*  In  InelasClc  Collisions 
1 .  H2(v)  +  >  H^(v  <  v)  +  Ny 

Due  Co  Che  face  ChaC  Che  race  for  Chls  process  Is  known  only  for  v 
»  1  ac  room  CemperaCure,  we  muse  cheoreClcally  exCend  Chese  resales  Co  Che 
lower  cemperacures  and  higher  levels  of  exclcaClon  encounCered  In  Che  upper 
acmosphere.  Simple  meChods  for  esclmacing  races  of  vlbraClonal-Co-vlbraClonal 
(V~V)  and  vlbraclonal  (V-T)  energy  cransfer,  such  as  SSH  cheory  and  Ics 
derlvaclves  (Lamberc  1977  and  references  chereln)  give  unreasonable  resales 
for  ^2  thus  may  not  be  used. 

Some  aspects  generally  seen  in  energy  transfer  collisions  are 
expected  be  true  for  this  process.  Specifically,  we  expect  all  such 
deacClvaClons  Co  occur  In  single  quantum  transitions  and  Che  deactivation  rate 
to  be  proportional  to  the  level  of  H2  excitation: 


V,  v-1 


V  k 


1,0 


(14) 


Rather  Chan  assuming  the  usual  T  dependence  for  deactivation  rates  (Lambert, 
1977)  in  this  system,  we  will  assume  that  the  k^  are  T  independent.  Our 
motivation  for  this  comes  from  noting  that  in  order  to  ensure  that  transition 
probabilities  obey  microscopic  reversibility,  simple  expressions  for  them  must 
contain  a  term  of  the  form  exp(AE/2kT).  AE  is  the  energy  being  converted  from 
vibration  to  translation  and  is  positive  for  deactivation  and  negative  for 
excitation.  Combining  this  dependence  with  usual  T  dependence  we  see  in  the 
simple  theories 


oc  exp(-BT  )exp(AE/2kT) 


(15) 
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Table  2:  State  Specific  Rates  k^(T)  for  the  Reaction  0  +  H2(v)  -*■  OH  +  H 


V 

a'*'  b  (°K) 

V  V 

k^(200K) 

k^(300K) 

0 

1.72(-9) 

740 

5220 

5.85(-22) 

3.28(-18) 

.  1 

2.32(-9) 

703 

2850 

1.11(-16) 

1.21(-14) 

2 

2.78(-9) 

650 

1820 

2.35(-14) 

4.72(-13) 

.  .  3 

3.37(-9) 

556 

1210 

6.48(-13) 

4.90(-12) 

4 

3.69(-9) 

512 

1010 

1.96(-12) 

l.ll(-ll) 

5 

4.04(-9) 

460 

800 

6.43(-12) 

2.63(-ll) 

*Rate 

constants  have 

units  of 

3  -1-1 

cm  molec  sec  . 

The  nimibers 

in  parentheses 

represent  the  power  of  ten  by  which  the  non-enclosed  number  should  be 
multiplied. 
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Clearly,  at  sufficiently  low  T,  Pj^q  will  increase  with  T.  Using  reasonable 
values  for  the  collection  of  parameters  B  and  AE  we  expect  no  enormous 
reduction  in  Pj^q  from  300-200K  (no  more  than  a  factor  of  three),  and  thus 
assume  no  T  dependence  for  the  deactivation  rate.  This  assumption  probably 
leads  to  a  slight  overestimate  of  the  ^2  quenching  rate  and  thus  a  slight 

•k 

underestimation  of  concentrations  and  hence  of  its  role  in  atmsopheric 
chemistry. 

Since  our  neglect  of  the  transmission  coefficient  in  section  B  above 
probably  led  to  an  overestimate  of  the  rate  of  R2b,  we  think  these  assumptions 
should  more  or  less  cancel  each  other  out. 

2.  H»)(v)  -I-  0  ->•  H2(v"<  v)  -h  Q 

It  is  known  that  vibrational  deactivation  in  potentially  reactive 
collisions  may  be  an  extremely  rapid  process  (Smith,  1976),  and  for  this 

A 

reason  we  must  consider  this  process  as  a  potential  loss  mechanism  for  H2  . 
In  addition,  it  has  been  seen  in  collinear  quantum  mechanical  scattering 
calculations  in  atom-dlatomlc  molecule  systems  that  multi-quantum  transitions, 
may  account  for  a  substantial  fraction  of  deactivation  collisions  (Kaye, 
1982). 

In  order  to  assess  the  relative  magnitude  of  this  process  and,  in 
particular,  the  Importance  of  multi-quantum  transitions  we  have  performed 
collinear  quantum  mechanical  scattering  calculations  on  this  system  using  the 
MOD  POL-CI  RMOS  surface  (Lee  et  al. ,  1982a).  For  simplicity,  we  will  consider 
one  n-quantum  transition  as  n  single-quantum  transitions.  We  estimate  the 
rate  of  deactivation  by  this  process  to  be  given  by 

k  (T)  -  R  (T)k^^^(T)  (16) 

v,v-l  V  V 

(2) 

where  the  superscript  represents  the  reaction  number;  k^  is  given  by  eq.  12 
or  alternatively  by  eq.  13  and  Table  2,  and  Ry(T)  is  determined  from  the 
collinear  quantum  mechanical  calculations  by  the  expression 

R^(T)  (v-v')k^^^  (T)  (17) 

V  »  1 _ vv _ 

k^^^(T) 

V 

where  the  k's  in  eq.  17  are  the  calculated  collinear  ones.  Results  of  the 
calculations  indicate  nearly  T  independent  values  of  Ry(T)  of  0.4,  0.7,  and 
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0.9  for  V  »  1,2,  and  3  respectively.  We  further  assume  values  of  1.0  for  v  « 
4  and  5.  The  assumption  that  one  may  use  the  results  of  collinear 
calculations  to  make  predictions  about  the  real  three-dimensional  world  has 
been  made  frequently  in  the  past  (Bowman  £t_  al. ,  1981;  Lee  al.  ,1982b; 

Wassam  and  Lee,  1976;  Wassam,  1980).  and,  while  certainly  not  being  exact, 
should  give  physically  reasonable  results. 

IV.  Aeronomy  of 

To  evaluate  the  possible  importance  of  R2b  in  controlling  the  H2 
concentration  in  the  earth's  upper  atmosphere,  we  will  consider  a  simple 
steady-state  purely  photochemical  model  (i.e.  containing  no  transport)  for  the 
80  km  altitude  region,  at  which  H2  production  is  expected  to  be  a  maximum 
(Hunten  and  Strobel,  1974;  Liu  and  Donahue,  1974).  For  convenience  we  adopt 
results  from  the  former.  Although  numerous  reaction  rate  constants  and 
branching  ratios  have  been  assigned  new  values  since  that  work  was  performed, 
the  qualitative  picture  remains  unchanged.  We  thus  assiime  that  the 

ii 

conclusions  we  draw  about  the  role  of  should  also  be  valid  in  a  model  using 
more  recent  values  for  the  many  parameters  necessary. 

Concentrations  of  important  species  and  reaction  rates  are  taken  from  the 
appendix  and  Figures  2  and  3  of  Hunten  and  Strobel  (1974)  and  are  shown  in 
Table  3  here.  Pertinent  rate  constants  and  reaction  rates  are  shown  for  v  » 
0-5  at  80  km  in  Table  4. 

* 

The  assumption  of  steady  for  H2  implies  the  following  expressions: 


-  0  -  R  +  R<’>,r0i  + 

dt  V  V  ^  V  V  <■  - 


k^^^^  FmI  +  k^^^[0(lD)l 1!H2(v) 


+(k^J[M]  +  k^J^  [01)[H2(v+1)1  (18) 

where  we  repeat  that  we  are  considering  only  single  quantum  deactlvatlons  and 
we  neglect  vibrational  excitation.  These  equations  may  be  rewritten  in  matrix 
form 


Kx  - 


R 


(19) 
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* 

Table  3;  Parameters  of  Importance  to  Aeronomy 


(Taken  from 

appendix  and  Figures  2 

and  3  of  Hunten  and 

Strobe 1. 

1974) 

z(km)  T(K) 

[^1 

[0^ 

f0(lD)l 

[H2I 

,(3a)^ 

50 

271 

2.14(16) 

3.25(9) 

3.05(2) 

1.5(9) 

1 

0 

60 

243.3 

6.47(15) 

3.25(9) 

1.3(2) 

2(9) 

2 

0 

70 

216.6 

1.70(15) 

2.0(9) 

4.8(1) 

1.5(9) 

6 

25 

80 

186 

3.76(14) 

7.0(9) 

3.5(1) 

8(8) 

80 

100 

90 

183 

6.39(13) 

2.2(11) 

1.5(2) 

1.5(8) 

0.05 

6 

100 

193.5 

9.74(12) 

8.0(11) 

6.0(2) 

2(7) 

0.001 

0.1 

*  Concentrations  are  in  units  of  cm  ,  Reaction  rates  (R)  are  in  units  of 

—  1  —1 

cm  sec  .  Numbers  enclosed  in  parentheses  are  the  powers  of  10  by  which 
the  non-enclosed  numbers  should  be  multiplied. 


Table  4;  Data  for  Steady  State  Calculation  at  80  k.m* 


V 

R(4a) 

k<2) 

k(^) 

r(7> 

r('> 

r(7> 

r(5) 

V 

V 

V 

V 

V 

V 

V 

0 

28 

/.  /. 

8.4(-24) 

0 

0 

5.9(-l4) 

0 

6.7(-9) 

I 

20 

24 

3.9(-17) 

1.6(-17) 

l.l(-7) 

2.7(-7) 

2.4(-2) 

6.7(-9) 

2 

14 

14 

1.2(-14) 

8.4(-15) 

5.9(-5) 

8.4(-5) 

4.8(-2) 

6.7(-9) 

3 

10 

8 

4.2(-13) 

3.8(-13) 

2.6(-3) 

2.9(-3) 

7.2(-2) 

6.7(-9) 

4 

7 

5 

1.4(-12) 

1.4(-12) 

9.8(-3) 

9.8(-3) 

9.6(-2) 

6.7(-9) 

5 

0 

6 

4.8(-12) 

4.8(-12) 

3.4(-2) 

3.4(-2) 

^4 

I 

CM 

• 

6.7(-9) 

*  Reaction  rates  (R)  are  in  units  of  cm  sec  ,  rate  constants  are  in  units 
3  —1  —I 

of  cm  molec  sec  .  Numbers  enclosed  in  parentheses  are  the  powers  of  10 
by  which  the  non-enclosed  numbers  should  be  multiplied. 


14 


where  x  is  the  column  vector  of  H2(v)  concentrations,  R  is  the  column  vector 
of  production  rates,  and  K  is  the  transition  matrix  whose  elements  are 
given  by 


K  ' 
vv 


fOl  +  +  k^^^ro(iD)rL  v' 

k^_^?[Ml  +  k^^?  I’Ol 


‘v+li 


v+1 


V 

v  +  1 


(20) 


Equations  18-20  do  not  include  the  effects  of  transport,  which  will  only 
be  important  for  H2(''"^)*  This  neglect  implies  that  our  calculated  values  of 

the  concentration  of  H2(v*0)  will  be  too  high,  but  this  error  will  not  affect 

* 

our  calculated  H2  destruction  rate  by  processes  involving  H2  •  The  calculated 
values  for  rH2(v)1  and  the  loss  rates  due  to  reaction  with  0(^P)  are  shown  in 
Table  5. 

The  major  result  to  come  from  this  study  is  seen  at  the  bottom  of  Table 

■k 

5:  the  total  loss  rate  of  H2  expected  due  to  chemical  reactions  of  H2  with 

0(^P)  is  approximately  3.9  cm”^sec'’^,  compared  with  a  total  H2  production  rate 
(see  Table  2)  of  180  cm“^sec“^.  Thus,  reactions  of  H*  remove  only  some  17,  of 
the  total  H2  produced  at  this  altitude  (80  km).  The  only  other  chemical  loss 
process  for  H2  important  at  80  km  is  R5,  which  occurs  with  a  rate  of 
5.4  cm^^sec”^.  Thus,  while  ™akes  a  substantial  contribution  to  chemical 
loss  of  H2,  the  total  contribution  of  chemical  losses  to  the  overall  H2  budget 
is  small;  most  of  the  H2  produced  at  80km  is  transported  to  higher  or  lower 

levels  of  the  atmosphere  (Hunten  and  Strobel,  1974).  The  very  small 

* 

contribution  of  loss  of  H2  by  reaction  of  H_  to  the  total  H2  budget  suggests 

^  (2) 

that  even  if  we  have  grossly  underestimated  k^  ,our  conclusion  will  still  be 
valid. 

From  Table  2,  one  may  see  that  80  km  is  the  most  favorable  altitude  for  a 
contribution  of  H*  to  the  overall  H2  chemistry.  Away  from  80  km,  H2 
production  is  reduced,  immediately  reducing  the  potential  role 
of  H*  chemistry.  Also,  [0(^0)]  is  a  minimum  at  80  km  and  thus  at  other 
altitudes  chemical  loss  of  H2  by  R5  will  be  greater.  At  lower  concentrations 
the  vibrational  quenching  rates  will  increase,  further  limiting  the  possible 
contribution  of  H*  .  At  higher  altitudes  the  combination  of  greatly  decreased 
H2  production  and  increased  [0(^0)]  will  approximately  cancel  the 
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★  1 

Table  5:  Contributions  of  H2  to  Loss  of  H2  Reaction  with  0(^P)  at  80  km* 


V 

,  (2),  3  ,  -1  -1. 

k  (cm  molec  sec  ) 

[H2(v);(cm  ^) 

V 

0 

8.4(-24) 

8(8)^ 

4.7(-5) 

L 

3.9(-17) 

^.3(3) 

1.2(-3) 

2 

1.2(-14) 

1.3(3) 

l.l(-l) 

3 

4.2(-13) 

4.3(2) 

1.3(0) 

4 

1.4(-12) 

1-5(2) 

1.5(0) 

5_ 

4.8(-12) 

3.2(1) 

1.1(0) 

sum 

— 

— 

3.9(0) 

i  Numbers  in  parentheses  are  the  powers  of  10  by  which  the  non-enclosed 
numbers  should  be  multiplied. 

^  ■a2(v  ■  0)1.  is  taken  from  Hunten  and  Strobel  (1974)  and  not  from  the 

results  of  our  calculation,  which  neglects  transport. 
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increased  "0(^P)^  and  decreased  total  pressure,  which  would  tend  to  increase 

*  * 

the  importance  of  reactions  of  •  Thus,  it  is  apparent  that  chemistry 

will  not  significantly  perturb  the  H2  budget  at  any  altitude  in  the  mesosphere 

or  lower  thermosphere.  We  emphasize  the  difference  between  our  results  and 

those  of  the  more  complete  but  less  chemically  accurate  model  of  Crutzen 

(1974),  who  suggested  that  consideration  of  R2b  might  greatly  reduce  upper 

atmospheric  H2.  Our  differences  arise  mainly  from  the  different 

(2) 

approximations  made  for  the  effect  of  vibrational  enhancement  on  .  The 

value  assumed  by  Crutzen  (1976)  for  v  >  1  is  greater  than  our  estimated  value 

for  V  <  5  at  80  km.  Since  our  model  contains  a  much  more  realistic  chemical 

(2) 

treatment  for  the  effect  of  H2  vibrational  excitation  on  k^  ,  we  feel 
confident  that  our  results  more  accurately  represent  the  effect 

A 

of  H2  chemistry  on  atmospheric  H2» 

V.  Conclusions 

We  have  shown  using  a  realistic  chemical  model  that  the  contribution  of 

chemical  reactions  of  vlbrationally  excited  H2  to  the  total  upper  atmospheric 

H2  budget  is  negligible.  This  conclusion  differs  greatly  from  the  preliminary 

result  of  Crutzen  (1974)  who  suggested  such  a  process  could  lead  to  a  dramatic 

reduction  in  upper  atmospheric  H2.  Our  model  contains  chemically  reasonable 

values  for  state-specific  reaction  rates  and  branching  ratios  for  the 

production,  chemical  reaction,  and  quenching  of  H.  .  Our  estimate  that 

^  * 

only  '.2%  of  the  H2  produced  at  80  km  is  lost  by  reaction  of  H2  is  sufficiently 
small  that  the  conclusion  should  be  independent  of  the  exact  assumptions  made 
in  our  model.  Gross  errors  (l.e.  more  than  an  order  of  magnitude  would  be 
necessary  to  qualitatively  alter  this  conclusion. 
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